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This is a review of current theory of black-hole dynamics, concentrating on the framework in 
terms of trapping horizons. Summaries are given of the history, the classical theory of black holes, 
the defining ideas of dynamical black holes, the basic laws, conservation laws for energy and angular 
momentum, other physical quantities and the limit of local equilibrium. Some new material concerns 
how processes such as black-hole evaporation and coalescence might be described by a single trapping 
horizon which manifests temporally as separate horizons. 
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I. INTRODUCTION 

Black holes are now generally regarded as astrophysi- 
cal realities, which are expected to be major sources of 
gravitational radiation, prompting extensive studies of 
dynamical, strong-field processes such as binary merg- 
ers. The textbook theory of black holes, however, mostly 
concerns stationary black holes or physically unlocatable 
event horizons^ 2 -' 3 -'^. In recent years, a new paradigm 
for dynamical black holes has been developed in terms 
of trapping horizons, generated from surfaces where light 
is momentarily caught by the gravitational field. They 
locate the black hole in a practical, local way. 

The original idea due to Penrose 5 - is that of a trapped 
surface, which compose the interior of a Schwarzschild 
black hole. Such surfaces played a key role in the singu- 
larity theorems of Penrose and Hawking 5 -^. The bound- 
ary of the black hole is composed of marginally trapped 
surfaces, here called marginal surfaces. They are also 
commonly called apparent horizons, though the textbook 
definition of the latter differs^. By itself, this does not 
capture the idea of a black hole, as such surfaces also exist 
in cosmological models. Presumably for this reason, such 
quasi-local ideas were not much developed at the time. 
Instead history took a different turn, with the introduc- 
tion of global methods and the near-universal acceptance 
of an event horizon as the definition of black hole. 

The author's interest in this area began with a rather 
simple observation: that it is outgoing (rather than in- 
going) wavefronts which are marginally trapped, that in- 
going wavefronts are converging, and that the trapped 
surfaces are inside (rather than outside). This led to a 
definition of a black hole by a future (rather than past) 
outer (rather than inner) trapping horizon 8 ' 9 . Here trap- 
ping horizon is used to mean any hypersurface foliated by 
marginal surfaces, a change from the original definition. 

Remarkably, this refinement sufficed to derive some 
key expected properties of black holes, assuming positive- 
energy conditions: the marginal surfaces have spherical 
topology; the horizon is spatial or null, and null only 
when a certain energy density vanishes, so that the hori- 
zon is one-way traversable; and the area A is constant 
when null and increasing when spatial. Physically, if 
something falls into a black hole, the horizon moves out 



and the area increases. 

In the case of spherical symmetry, a rather complete 
picture emerge d 1 ^ 11 ' 12 . There is a standard definition 13 
of active gravitational mass M. It satisfies an energy 
conservation law akin to the first law of thermodynam- 
ics, with energy-supply and work terms, reducing to the 
Bondi energy equation at null infinity. There is also a 
natural definitioni 1 - of surface gravity k, taking a quasi- 
Newtonian form, just M/r 2 in vacuo, in units G = 1, 
where A = Airr 2 . It has various desirable properties, in- 
cluding that k > 0, k — or k < on outer, degenerate 
or inner trapping horizons respectively. Then projecting 
the energy conservation law along a trapping horizon, de- 
noted by a prime, yields M' = kA' /8ir plus work term, 
where = denotes evaluation on a trapping horizon. This 
is a dynamic version of the so-called first law of black- hole 
mechanics^. 

Versions of these definitions and results also exist in 
cylindrical symmetry^, which is the simplest situation 
allowing gravitational radiation, in plane symmetry 15 
and in a quasi-spherical approximatio n 16 i 17 ' 18 i 19 , which 
allows the interaction of roughly spherical black holes 
with gravitational radiation. In all these cases, there is 
an effective energy tensor Q a /3 for the gravitational ra- 
diation, entering the energy conservation law additively 
with the matter energy tensor T a p. 

Meanwhile, Ashtekar and others developed a theory 
of isolated horizons^^^ 2 ^* 2 ^^, which are null 
trapping horizons with a hierarchy of additional condi- 
tions, each intended to describe a black hole in local equi- 
librium to some degree, generalizing Killing horizons. In 
particular, they defined angular momentum J, entering 
a first law in the expected way. 

A breakthrough came with the work of Ashtekar & 
Krishnan 29 ,30 ! 31 on what they called dynamical horizons, 
which are spatial future trapping horizons. They derived 
flux laws for both energy and angular momentum. In the 
versions subsequently developed by the autho r 32 ' 33 i 34 ' 35 , 
they hold for any trapping horizon and take the form of 
conservation laws 

L e M S j> s *{T aP + @ a0 )k a r p (1) 
L e J - j> s *{T a p + Q a p)^ a T P (2) 
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where S is a marginal surface, *1 its area form, L is 
the Lie derivative, £ is a generating vector of the trap- 
ping horizon, r its normal dua l 36 ' 37 ' 38 , and k and -0 are 
certain vectors which play the role normally played by 
Killing vectors. Here Q al j is again interpreted as an effec- 
tive energy tensor for gravitational radiation. Apart from 
the inclusion of a /3, these have the same form as corre- 
sponding flat-space conservation laws in surface-integral 
form. 

Some subtleties remained. Firstly, in what one might 
call the bottom-up approach of Ashtekar and co-authors, 
the formalisms for isolated and dynamical horizons were 
quite different, and it was not clear how to relate them. 
Some progress was made by Booth & Fairhurst for slowly 
evolving horizon o 39 ' 40 ^ 41 . Conversely, in the author's top- 
down approach, whereby local equilibrium is given simply 
by the trapping horizon being null, there is a degeneracy 
in taking the null limit. Remarkably, consideration of an- 
gular momentum closed this gap by suggesting a natural 
way to fix the degenerac y 34 i 35 . The resulting gauge-fixed 
null trapping horizon is then equivalent, modulo further 
gauge-fixing, to a weakly isolated horizo n 24 i 25 i 26 i 27 . Thus 
the top-down and bottom-up approaches have converged 
on the same definition of a black hole in local equilibrium. 

For this reason, combined with the existence of the 
reviews of Ashtekar & Krishnan^i, Boothi 2 -, Krishnan^ 3 - 
and Gourgoulhon & Jaramillo 4 ^, this review will con- 
centrate on the paradigm in terms of trapping horizons, 
without attempting to make detailed comparisons with 
the other approaches. Indeed, it is intended to be largely 
complementary to those reviews, apart from inevitable 
overlap on core issues. No attempt will be made at a 
comprehensive list of references, since the area and re- 
lated areas are rapidly evolving. With a broad audience 
in mind, detailed calculations and proofs will be deferred 
to references. General Relativity, including the Einstein 
equation, will be assumed throughout, though the essen- 
tial ideas can be generalized to other metric theories or 
dimensions. 



II. CLASSICAL THEORY OF BLACK HOLES 

What one might nowadays call Newtonian black holes 
were described long ago by Mitchell and Laplace 2 -, as 
objects for which the escape speed \j2GMjr is greater 
than the speed c of light. Any light leaving such an object 
would be trapped by the gravitational field and fall back, 
rendering it invisible. Curiously, the Newtonian relation 
r < 2GM/c 2 will survive quite generally in a sense to be 
revealed, with units c = 1 = G henceforth. 

Almost immediately after Einstein formulated General 
Relativity^ 5 - in 1915, Schwarzschild 46 derived the solution 
for the gravitational field of a point mass M. Charge Q 
was soon added by Reissner— and Nordstrom 48 . How- 
ever, the geometry of the solutions was not understood 
for decades. Einstein & Rosen~ deciphered the (Q = 0) 
spatial geometry in 1935, with a wormhole connect- 
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FIG. 1: Example of the unknowable nature of event horizons: 
Penrose diagram of a spherically symmetric space-time which 
is initially flat but contains an ingoing radiation shell, forming 
a Schwarzschild black hole. An observer in the flat region 
crosses the event horizon but feels nothing. 



ing two asymptotically flat spaces, but it was not until 
around 1960 that the maximally extended space-time ge- 
ometry was understood, with work of Wheeler, Kruskal 50 
and Fronsdal 5 - 1 -. The space-time diagram clearly shows a 
region r < 2M from which light cannot escape. It had 
already been shown by Oppenheimer & Snyder— in 1939 
that gravitational collapse could indeed produce such a 
region. According to relativistic causality, nothing else 
can escape either. 

Angular momentum J was added by Kerr— and (for 
Q 0) Newman et ali^ 4 - and finally around 1968, the 
term "black hole" entered the lexicon, attributed to 
Wheeler—. A few years of rapid progress followed, culmi- 
nating in the classical paradigm. Black holes were defined 
in general by event horizons as defined by Penrose^, 
which satisfy the area-increase theorem of Hawking 57 , 
A' > 0. This became known as a second law, part of the 
four laws of black-hole mechanics formulated by Bardeen, 
Carter & Hawkingi, which seemed analogous to the four 
laws of thermodynamics. The zeroth law is that surface 
gravity k is constant on stationary black holes. The first 
law is 

5E = kSA/8tt + QSJ + $6Q (3) 

for perturbations of stationary black holes, where E is the 
ADM energy at spatial infinity, 17 is the angular speed 
and $ is the electric potential. The third law is n 0, 
by perturbations of stationary black holes. The discov- 
ery of quantum black-hole radiance by Hawking^ 8 -, with 
temperature n/2ir in units k = 1, led to much research 
on black holes as a key area to generate and test ideas 
concerning the interface of gravity, quantum theory and 
thermodynamics. 

Once a theoretical curiosity, observational evidence for 
black holes has by now accumulated past the point where 



3 



they are part of mainstream astrophysics, both as stellar- 
mass black holes produced as supernova remnants, and 
as supermassive black holes mysteriously found at the 
cores of most galaxies, powering active galactic nuclei. 
Gravitational-wave detectors are expected to be able to 
observe such processes as binary inspiral and merger of 
black holes or neutron stars, prompting much research in 
numerically computing waveforms. 

In such dynamical situations, the classical paradigm 
is of little help. Firstly, event horizons are defined as 
the boundary between where light can or cannot escape 
to infinity, which means not only infinite distance but 
infinitely far in the future. Therefore one cannot locate 
it without knowing future events, indeed the entire future 
history of the universe, which is beyond mere mortals. It 
is easy to give examples*^ where an event horizon exists 
in a region of flat space-time (Fig.[T]). Then an observer 
crossing the event horizon cannot detect it in any way 
and has never experienced a gravitational field. 

In practice, it is normally much easier to locate 
marginal surfaces in numerical simulations. Note also 
that three of the above laws are more accurately de- 
scribed as black-hole statics, since they concern station- 
ary black holes, specifically Killing horizons. So the ques- 
tion arises: are there corresponding laws of black- hole 
dynamics? 



III. DYNAMICAL BLACK HOLES 

The basic idea is easily explained. Imagine enclosing 
a star with a roughly spherical spatial surface at some 
moment of time, and detonating a flash of light simul- 
taneously at each point of the surface. Two wavefronts 
form, one ingoing and one outgoing. Normally one ex- 
pects the outgoing wavefront to have increasing area and 
the ingoing wavefront to have decreasing area. This is 
measured at each point by the expansions 9± of light rays 
in the wavefronts: 9+ > for the outgoing wavefront and 
< for the ingoing wavefront. However, the gravita- 
tional field of the star tends to drag things toward it, 
including light. Thus the outgoing wavefront does not 
expand as much as if the star were not present. The ef- 
fect increases closer to the star and for larger mass, as 
well known from Newton's law of gravitation. For large 
enough mass, for a surface close enough, it may happen 
that the outgoing wavefront has decreasing area. This is 
a black hole: all light and therefore matter are confined 
inside a shrinking area. 

The idea then is that the boundary of the black hole 
at a given time is a marginal surface, meaning that one 
wavefront has instantaneously parallel light rays, in this 
case 9 + — 0. To characterize a black hole, it is also 
important that the ingoing wavefront is converging, < 
0, and that 9+ is decreasing in the ingoing direction. 

This is formulated geometrically as follows. A spatial 
surface 5* in space-time has two unique future-pointing 



null normal directions, along null normal vectors l±: 

g(i±,i±) = Q, n ± = o (4) 

where g is the space-time metric and _L denotes projec- 
tion onto S. The null expansions are 

9±=L ± \og*l (5) 

where L± = L\ ± . 

Then S is said to be trapped if 9 + 9- > 0, marginal if 
9 + = or L =0, and untrapped (or mean convex) if 
9+9- < 0. In terms of the expansion vector or mean- 
curvature vector 

H = g- 1 (dlog*l) (6) 

this is equivalent to H being temporal, null or spatial, 
respectively. One finds 

H = -e/(9-l + + 6+lJ) (7) 

where / is a normalization function, 

e-f = -g(l+,l-). (8) 

For any normal vector 77, 1.7] = 0, the expansion 9 V = 
L n log *1 is given by 8 V — g(H, rj). 

Untrapped surfaces have a local spatial orientation: an 
achronal (meaning spatial or null) normal vector 77 is out- 
ward or inward if 9 V > or 9 V < respectively. Locally 
one can conventionally fix 9+ > 0, < in an un- 
trapped region; then l + is outward and Z_ is inward. Con- 
versely, trapped surfaces have a local causal orientation: 
if H is future or past causal (meaning temporal or null), 
the surface is future or past trapped respectively. Then 
future or past trapped surfaces have 9± < or 6± > 
respectively. Marginal surfaces have both orientations, if 
the other null expansion has fixed non-zero sign. 

A trapping horizon is a hypersurface foliated by 
marginal surfaces. It is said to be outer or inner if 
L-9 + < or L-9 + > respectively, for the case 9 + = 0, 
where the null normals l± are extended off the horizon 
to generate ingoing and outgoing wavefronts from the 
marginal surfaces, i.e. two families of null hypersurfaces 
labelled by , intersecting in the marginal surfaces, such 
that lA(dx B ) = 6 a- This is locally unique unless the 
trapping horizon is itself null. Equivalcntly, the horizon 
is outer or inner if V- (H — H*) > or V ■ (H - H*) < 
respectively, where rj* = rj + l+ — rj~l- is the normal dual 
vecto r 36 i 37 i 38 to a normal vector rj = i] + l+ + r)~l-: 

_L?y*=0, g(rj*,Tj) = 0, g(rj* ,rj*) = -g(rj,Tj). (9) 

This is a vector version of the normal Hodge dual for 
1-forms. 

Then a future (respectively past) outer trapping hori- 
zon provides a local definition of a generic black (respec- 
tively white) hole£. More precisely, the idea is that a 
non-degenerate black hole exists only if such a horizon 
exists. As to the converse, see the Concluding Remarks. 
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Note that outer versus inner has been defined with re- 
spect to the ingoing null direction, as this is invariant, 
rather than in some ingoing spatial direction. Also, de- 
manding a strict sign for L-8+ has excluded extremal 
black holes, for which L_6* + vanishes. This is for simplic- 
ity only, as they can be treated as special cases. Relaxing 
the definition to L-6+ < would allow cases which are 
not black holes. 



IV. BASIC LAWS 

Summarized here are some basic properties of trapping 
horizons 8 . Where not obvious, the proofs involve either 
(for 8+ = 0) the T ++ component of the Einstein equa- 
tion, where the null energy condition (NEC) is cited, or 
the component, where the dominant energy condi- 
tion (DEC) is cited. One introduces a normal generating 
vector £ of the marginal surfaces in the horizon, which 
therefore satisfies L^9+ = on the horizon. The area of 
the marginal surfaces, if compact, is 

A = <j> *1. (10) 

Trapping: for a future/past, outer/inner trapping hori- 
zon, there are trapped surfaces to one side and untrapped 
surfaces to the other side. For a future outer trapping 
horizon, this reflects the defining idea that outgoing light 
rays are momentarily parallel, 6 + = 0, diverging just out- 
side, 6 + > 0, and converging just inside, 6+ < 0, while 
ingoing light rays are converging, < 0. 

Signature: assuming NEC, an outer or inner trapping 
horizon is achronal or causal respectively, and null if and 
only if the effective ingoing energy density 7++ + 6++ 
vanishes, where the explicit expressions for Q a /3 are given 
later (|21j) . In particular, this means that black-hole hori- 
zons are one-way traversable: one can fall into a black 
hole but not escape, at least through the outer horizon. 

Area: assuming NEC, future outer or past inner trap- 
ping horizons have non-decreasing area form, 8^ > 0, and 
therefore (if compact) non-decreasing area, L^A > 0, in- 
stantaneously constant (8$ = 0) if and only if the horizon 
is null. For past outer or future inner trapping horizons, 
all signs reverse. Here the orientation of £ is such that, 
in the null limit, it is future-null, which means that it is 
future causal for inner horizons and outward achronal for 
outer horizons. In particular, this means that black holes 
grow if they absorb any matter or gravitational radiation, 
and otherwise remain the same size. 

Topology: assuming DEC, a future/past outer trapping 
horizon has marginal surfaces of spherical topology (if 
compact). The proof uses the Gauss-Bonnet and Gauss 
divergence theorems. Thus realistic black holes are topo- 
logical^ spherical. If degenerate horizons are considered, 
then toroidal topology is just allowed, but highly non- 
generic, in particular Gaussian flat, and so presumably 
unstable. 



Area limit: assuming DEC and a positive cosmological 
constant A, outer trapping horizons satisfy^ A < Att/A. 
Thus black holes are smaller than the cosmological hori- 
zon scale, corresponding to an area 127r/A. 

Pedagogically, this would be the place to give 
details in spherical symmetr y 10 ! 11 ' 12 , cylindrical 
symmetry^, plane symmetry^ and the quasi-spherical 
approximatio n 16 ' 17 ' 18 ' 19 , but they are omitted here due 
to limitations of space, apart from as explained in the 
Introduction. 



V. CONSERVATION OF ENERGY 

The simplest generalization of the Schwarzschild rela- 
tion 1 — 2M/r = g rr to a general surface is the Hawking 
mass 6 -^ 

M= r -(l-±-j s , 9{ H,H)) (11) 

= r 2{ i+ h£* efe+e -) 

where 

r = ^fA/An (12) 

is the area radius. It has various useful properties. Large 
spheres: in an asymptotically flat space-time, M tends 
to the Bondi or ADM energy at null or spatial infinity, 
respectively. Small spheres: M/volume — > density at a 
regular centre. Trapping: a surface is trapped, marginal 
or untrapped if r < 2M, r = 2M or r > 2M respectively. 

In particular, this means that M is the irreducible mass 
of a future outer trapping horizon, L^M > 0, assum- 
ing NEC. This follows directly from the area law, since 
A = 4ir(2M) 2 , where = henceforth denotes evaluation 
on a marginal surface. Recall the irreducible mass for 
stationary black holesS: M = (|m(m+ (to 2 - a 2 ) 1/2 ) 1/2 
for Kerr black holes, the mass which remains even if ro- 
tational energy is removed by the Penrose process. The 
original concept arose due to quasi-stationary arguments, 
but here it is exact and gives a physical meaning to M. 

The simplest generalization of the Schwarzschild sta- 
tionary Killing vector k = d t is the canonical time 
vector— 

k = {g-^dr))* = e f {L + rl_ - L_rZ + ) (13) 

or equivalcntly 

±k = 0, k-dr = 0, g(k,k) = -g- 1 (dr,dr). (14) 

In particular, g(k,k) = and k = ±g^ 1 (dr) on a trap- 
ping horizon 8± = 0. So trapping horizons are char- 
acterized by k being null, just as Killing horizons are 
characterized by a stationary Killing vector being null 
(Figf2J). In spherical symmetry, k reduces to the Kodama 
vector—, whichi^ is a Noether current with Noether 
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FIG. 2: Penrose diagram of gravitational collapse to a black 
hole, satisfying cosmic censorship. When a trapping hori- 
zon (blue bold line) forms, the flow lines of k (magenta long 
dashes) and g~ 1 (dr) (green short dashes) switch over, both 
being null at the horizon. 




FIG. 3: A non-null hypersurface H foliated by spatial surfaces 
S, with generating vector £ and its normal dual r = £*. 

charge M, and is relatec—i to k by k^Vipk^ = ±Kk a 
on a trapping horizon 9± = 0, analogously to the usual 
definition of surface gravity for stationary black holes. 

Lastly, one needs the dual normal vector to the trap- 
ping horizon: 

r = r =£ + l+-Cl- (15) 

which can be locally chosen to be future-pointing for 
outward-pointing £. For a spatial trapping horizon, £ 
is spatial and r is temporal, while r — > £ as a trapping 
horizon becomes null (Fig|3J|. 

With these ingredients, conservation of energy takes 
the surface-integral for m 32 ' 33 

L € M^ £*(T AB + e AB )k A T B (16) 

where the components AB are given below. The proof is 

a calculation using the T++ and components of the 

Einstein equation and the Gauss-Bonnet and Gauss di- 
vergence theorems. It may be noted that evaluation on a 
trapping horizon yields remarkable cancellations in what 
would otherwise be lengthy expressions on the right-hand 
side. The conservation law can equivalently be written 
in volume-integral form 

[M] £* / *{T AB + e AB )k A T B A dx (17) 

where x labels the marginal surfaces, £ = d x , which ex- 
presses the change [M] in M between two marginal sur- 
faces in the horizon H. For a spatial trapping horizon 



with unit normal f — r / \jg xx and proper volume ele- 
ment il = *y/g X x A dx, there is the proper-volume form 

[M] = / i(T AB + e AB )k A f B (18) 

which has the same form as the usual expression for en- 
ergy, with k replacing a stationary Killing vector. How- 
ever f is ill defined and il — * if the trapping horizon 
becomes null, g xx — > 0, which is the physically important 
limit where a growing black hole ceases to grow, so the 
above surface-integral form is preferred. 

To spell out the components Q AB , introduce the trans- 
verse metric h ab , i-e. the induced metric of S, the null 
shears 

CT±ab = hJh s b L ± h lS - 9±h ab (19) 

and the normal fundamental forms 

C±a = e^±/sV 7 4). (20) 

Note here that indices a, (3 . . . are general, A, B . . . nor- 
mal, and a, b . . . transverse. Then the symmetric bilin- 
ear forms a± are transverse, a± — _Lcr±, and traceless, 
h ab a± a b — 0, while (± are transverse 1-forms, £± = _L£±. 
Then the components of Q a p with respect to l± are 

e±± = |k ± || 2 /327r, e ±T -e^|C±| 2 /8^ (21) 

where |C| 2 = h ab C a Q b > 0, ||cr|| 2 = h ab h cd a ac cr bd > are 
transverse norms. Here the signs indicate that Q a p satis- 
fies DEC, so that gravitational radiation carries positive 
energy. 

The components may be interpreted as gravitational 
energy densities, by geodesic deviation of test particles^: 
0++ is the ingoing transverse mode, reducing to the 

Bondi energy density at past null infinity; is the 

outgoing transverse mode, reducing to the Bondi energy 
density at future null infinity; |_ is the ingoing longitu- 
dinal mode, with r 2 |_ — >■ at null infinity; 0^ is the 

outgoing longitudinal mode, with r 2 0_| — > at null in- 
finity. The 0±± components also recover expressions for 
energy density of gravitational radiation, with 0±^ van- 
ishing, in the high-frequency linearized approximation 3 , 
in cylindrical symmetryii, plane symmetr}^ and in the 
quasi-spherical approximatio n 16 ' 17 ! 18 ' 19 . 

Physically, the energy conservation law expresses the 
increase in irreducible black-hole mass M in terms of the 
energy densities of the infalling matter and gravitational 
radiation. Since A = 4ir(2M) 2 , it also describes how a 
black hole grows. 

VI. CONSERVATION OF ANGULAR 
MOMENTUM 

The standard definition of angular momentum for an 
axial Killing vector ip and at spatial infinity is the Komar 
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integral^ 

J M = - IF" / *«a/jVV (22) 

where £,43 is the binomial. For a general transverse vec- 
tor ip, .Lip = ip, it can be rewritten a o 34 i 35 

J bP] = ^ £ (23) 

where the twistSS 

w« = i^ o/ ,p_,Z + f (24) 

is a transverse 1-form, _Lw = w, measuring the non- 
intcgrability of the normal space. 

For the weak-field metric in spherical polar coordinates 
(t, r, 1?, y>), J[9p] recovers the standard definition of angu- 
lar momentum 3 -. Also the precessional angular velocity 
(O-f)f — gf2 of a gyroscope in the unit direction f, due to 
the Lense-Thirring effect, is directly related to the twist 
by to ~ f2 x f. Thus the twist indeed encodes the twisting 
around of space-time due to a rotating mass. 

The twist is an invariant of a non-null foliated hyper- 
surface H, so the twist expression for J[ip] is also an in- 
variant of H. It coincides with the 1-form used to define 
angular momentum for dynamical horizons by Ashtekar 
& Krishnan 2 ^, but not that used for isolated horizons by 
Ashtekar et a h 23 i 24 , which is ui + \Df, where D is the 
covariant derivative of h. They will give compatible J[ip], 
by the Gauss divergence theorem, if the axial vector has 
vanishing transverse divergence, 

D a ip a £* 0. (25) 

If there exist angular coordinates (1?, ip) on S, completing 
coordinates (x, tp) on H, such that ip = d v , then recall- 
ing that £ = d x and that coordinate vectors commute, 

m> (26) 

which was previously proposed as a natural way to prop- 
agate ip along £ by Gourgoulhon^i. There have been var- 
ious suggestions for specifying ip more uniquel y 34 ' 35 i 68 ' 69 . 
Then conservation of angular momentum takes the 

form 34£5 

L^J = - ^*(T aB + ® aB )ip a T B . (27) 

The proof is a calculation using the T a B components of 
the Einstein equation, and requires various cancellations 
due to evaluation on a trapping horizon and the condi- 
tions p5 ]) -([26 ll . Here 

e a ± = --^-h cd D d a ±ac (28) 

is the transverse-normal block of the effective energy ten- 
sor for gravitational radiation. Recalling the energy den- 
sities 0±± = ||<7±|| 2 /327r (l2"Tj) . indicating that transverse 



gravitational radiation is encoded in null shear a~±, it 
seems that differential gravitational radiation has angu- 
lar momentum density. So this describes how a black hole 
spins up or down, due to infall of co-rotating or counter- 
rotating matter or gravitational radiation. 

Thus conservation of energy and angular momentum 
take a similar form JT]) , @ . Both take the same form as 
standard expressions in flat space-time for a stationary 
Killing vector k and an axial Killing vector ip, except 
for the inclusion of gravitational radiation in Q a p . They 
are the independent conservation laws expected for an 
astrophysical black hole, which defines its own centre-of- 
mass frame and its own axis of rotation. 

Higher source multipoles have been defined for isolated 
horizons by Ashtekar et aliSi and for dynamical hori- 
zons by Schnetter et al^. A gauge-dependent definition 
of linear momentum has been proposed by Krishnan et 
alj22, which may be useful in studying the recoil or kick 
effect in asymmetric binary mergers. 

VII. QUASI-LOCAL CONSERVATION LAWS 

To compare with the classical paradigm, one needs to 
include charge Q, which is defined in terms of charge- 
current density j as 

[<?] = -/ *<?(j, r) A dx = - [ *g(j, f) (29) 
Jh Jh 

where the more usual second expression holds only for 
a spatial hypersurface H. As above, the surface-integral 
form is 

L& = - f *g(j,r). (30) 

The above conservation laws can be written in the same 
form 

L t M S* - / *g(j, t), L^J = - <f *g{j, r) (31) 
Js Js 

by identifying current vectors 

f = -k A (T AB + e AB ), f = MT aB + e aB ). (32) 

The standard physical interpretation of the conserved 
vectors is j = (energy density, energy flux) , j — (angular 
momentum density, angular stress), j = (charge den- 
sity, current density). For spatial £: § s *g(j,£) — 
power; f s *g(jj,£) = torque; § s *g(j,Q = current; 
- §s*9(ji T ) = energy gradient; - § s *g{j,r) = angular 
momentum gradient; — <f s *g{j, t) = charge gradient. 
Charge conservation in local differential form is 

y a r = 0. (33) 

However, for energy and angular momentum, one has 
only quasi-local conservation laws holding on a trapping 
horizon: 

/ *v Q r = / *v Q r = 0. (34) 

Js Js 
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This subtly confirms the view that energy and angular 
momentum in General Relativity cannot be localized 3 , 
but might be quasi-localized as surface integrals, as long 
ago argued by Penrose 7 - 3 .. The corresponding conserva- 
tion laws have indeed been obtained in surface-integral 
but not local form. 



Here the state-space perturbations in the classical 
law for Killing horizons 1 , or the versions for iso- 
lated horizonsSBiSLSSiSSiSlSfiiSSiSl, have been replaced by 
derivatives along the trapping horizon, thereby promot- 
ing it to a dynamical law. 



VIII. STATE SPACE 



IX. EQUILIBRIUM: NULL TRAPPING 
HORIZONS 



There are now three conserved quanti- 
ties (M,J,Q), forming a state space for 
dynamical black holes. Following various 

au ^ nor : ,20. 21. 22. 23. 24. 25. 26. 27.28. 29. 30. 31. 39.40. 41.42. 43 re _ 

lated quantities may then be defined by formulas 
satisfied by Kerr-Newman black holes, specifically those 
for the ADM energy 



E = 



V((2A/) 2 +Q 2 ) 2 + (2J) 2 



AM 



the surface gravity 



(2M) 4 - (2 J) 2 - Q i 



2(2M) V((2M) 2 + Q 2 ) 2 + (2J) 2 
the angular speed 

nr J 



MV((2M) 2 + Q 2 ) 2 + 
and the electric potential 

((2M) 2 + Q 2 )Q 



{2Jf 



2M X /((2M) 2 + Q 2 ) 2 + (2J) 2 



(35) 



(36) 



(37) 



(38) 



These formulas may not be so familiar, since the classical 
theory uses (E, J, Q) as parameters with M dependent, 
but are easily obtained by inverting standard formulas. 

In the dynamical context, E > M is not the ADM 
energy, but can be interpreted as the effective energy of 
the black hole. Expanding for J <C M 2 and Q <C M, 
E ss M to leading order, then to next order, 



E « M + §/Q 2 + \Q 2 jr 



(39) 



where J = IQ, defines the moment of inertia I = 
My/((2M) 2 + Q 2 ) 2 + (2 J) 2 ^ Er 2 . Thus E includes ir- 
reducible mass M, rotational kinetic energy ss \l£l 2 and 
electrostatic energy « \Q 2 /r, the latter being standard 
Newtonian expressions. 
The state-space formulas 



K = 8lT 



dE 



1 dE 



dA 4M dM ■ 



^ dE ^ dE fArsS 
na fl7' (40) 



then yield a dynamic versio n 34 i 35 of the so-called first law 
of black- hole mechanics ([3]) : 



LfE = —LcA + QLtJ + $L 5 Q. 

07T 



(41) 



When a growing black hole ceases to grow, the generi- 
cally spatial trapping horizon becomes null. The central 
result here is that, assuming DEC, 

9(j»~9(j»~9(j^) = (42) 

on a null trapping horizo n 34 ' 35 . Therefore the conserved 
quantities are actually preserved: 

LtM L e J = L S Q = 0. (43) 

This indicates that local equilibrium is attained when a 
trapping horizon becomes null. 



It follows that 



L f E ^ Lck ^ L e fl ^ Lt<S> ^ 0. 



(44) 



In particular, the surface gravity, which already satisfies 
Dk = by definition (|36|) . is constant where a trapping 
horizon becomes null. This is a quite general zeroth law. 
The above result is stronger still, since it expresses com- 
plete local equilibrium, not just thermal equilibrium. 

By the area law&, which includes L^A = => H null, 
this also shows that a black hole cannot change its angu- 
lar momentum or charge without increasing its area. 

On a null trapping horizon, one may take £ = r = 
but the other null vector is non-unique, leading to 
non-uniqueness in lo, used to define angular momen- 
tum here and for dynamical horizons by Ashtekar & 
Krishna n 29 ' 30 ' 31 . However, uj + \Df is unique, an in- 
trinsic normal fundamental form of a null hypersurface, 
therefore used to define angular momentum for isolated 
horizons by Ashtekar et alM^^^^^^l^. On 
the other hand, the extrinsic normal fundamental form 
uj - \Df is preserved 3 ^: DEC =*> L+{lj - \Df) S 0. 
The non-uniqueness is therefore naturally fixed by 

Df = 0. (45) 

Recalling the definition ([8|) of the normalization function 
/, this is a legitimate choice of gauge. Then all three 
normal fundamental forms coincide. Also L^ui = 0, so 
that L{ J = assuming only that iji is a coordinate vector. 
Comparing with energy, NEC => L^M = automatically 
on a null trapping horizon. 

Thus consideration of angular momentum resolves the 
ambiguity in taking the null limit. The general formalism 
for trapping horizons then applies in all cases, describing 
transitions between growing and non-growing phases of 
a black hole. This has largely recovered the notion of 
weakly isolated horizon due to Ashtekar et al . 24 ' 25 ' 26 ' 27 , 
except that the (allowable) scaling freedom in £ has not 
been fixed 3 ^ 5 .. 
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FIG. 4: Black-hole formation and evaporation. The vertical 
coordinate is loosely an advanced time, while the horizontal 
coordinates loosely indicate area, with marginal surfaces in 
the outer and inner trapping horizons depicted respectively as 
blue and red circles. In the lower half of the picture, the black 
hole is forming and the outer horizon is growing, while in the 
upper half, the black hole is evaporating and the outer horizon 
is shrinking. The inner horizon has the opposite properties. 
No attempt has been made to indicate timescales. 




FIG. 5: Black-hole coalescence, with conventions as in Fig(4] 
No attempt has been made to indicate relative motion. 



X. TYPE-CHANGING HORIZONS 

While the above properties indicate that a future outer 
trapping horizon serves as a practical definition of black 
hole, it should be noted that a trapping horizon may 
change its type under evolution. For instance, in grav- 
itational collapse, it is common for an inner horizon to 
form simultaneously with an outer horizon. In fact they 
join smoothly and are really two parts of a single horizon 
in a space-time sense, manifesting as distinct horizons in 
a spatial slicing (Fig0]) . According to the basic laws, the 
transition occurs when the horizon is null, with the outer 
horizon achronal and the inner horizon causal. 

Booth et al^i constructed examples where an outer 
horizon becomes past-null, turning into a past-causal in- 
ner horizon, which can then turn into an outer horizon 
again, and so on. Also, numerical examples of Schnetter 
et ali^i show that a trapping horizon can be partly outer 
and partly inner on a given marginal surface, interpolat- 
ing between regions where it is strictly inner or strictly 
outer. 

One application of such ideas is to evaporating black 
holes. Ingoing Hawking radiation tends to have nega- 
tive energy density, T ++ < 0, violating NEC, and if it 
dominates over gravitational radiation in the sense that 
T ++ + ©++ < 0, the basic laws reverse: the outer hori- 
zon becomes causal and shrinks, while the inner horizon 
becomes spatial and grows. Thus it is possible that they 
simply reunite smoothly, closing off the region of trapped 
surfaces (FigHJ). The black hole has then evaporated. If 
no singularity ever formed, the space-time would have 
the global structure of Minkowski space-time and there 
would be no possibility of information loss and the as- 
sociated information puzzle. Explicit examples of such 



space-times have been constructed 7 -^ 

Another application is to the coalescence of binary 
black holes, of great interest as generators of strong 
gravitational waves. Typically a common outer hori- 
zon forms around the original two outer horizons. As 
above, it joins smoothly to an inner horizon, though the 
latter is not always tracked in simulations. It can be 
conjectured that the inner horizon later joins with the 
original outer horizons. For instance, the inner horizon 
may pinch off and change topology from one sphere to 
two spheres, which close in on the original outer hori- 
zons and merge smoothly with them, leaving the en- 
tire region inside the common outer horizon composed 
of trapped surfaces (Fig[5]). This appears to be consis- 
tent with simulations 71 , where the original outer hori- 
zons are slowly evolving and appear to be stable, while 
the inner horizon appears to be unstable and is rapidly 
evolving. More complex topologies are possible, and in- 
deed there is numerical evidence that the original outer 
horizons can osculate and then intersect each other—, al- 
lowing a self-intersecting inner horizon. It is also possible 
that singularities may intervene first, but this would be 
occurring in untrapped regions and so violate a version 
of strong cosmic censorship. 

According to the above picture, the entire set of hori- 
zons actually form a single trapping horizon in space- 
time, despite manifesting temporally as different and ap- 
parently unrelated horizons. The whole horizon can be 
smooth, though the spatial sections will have at least 
one point of non-smoothness where the topology changes. 
The irreducible mass M, being given by the area A, may 
be evaluated throughout the entire process, satisfying 
the energy conservation law. The separate areas A\ and 
A2 add to the common area A at the topology change, 
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A = A\ + A2, so that, assuming NEC, 

A > Ai + A 2 (46) 

at any stage where either exists, where one follows the 
original outer horizons outwards in a space-time sense 
until the trapping horizon becomes a common outer hori- 
zon. The inequality holds also during inner stages, where 
one should note that the inner horizon is past-causal 
in an outward sense, so the area is still increasing out- 
wards. This is an area-increase law for coalescing black 
holes, reminiscent of Hawking's famous theorem for event 
horizons^. 
It also implies 

M > sj Mf + Mf . (47) 

However, note that the total irreducible mass is ac- 
tually reducible in such a process, indeed it necessar- 
ily decreases discontinuously at the change of topology, 
M < Mi + Mi. The greatest fractional loss of total irre- 
ducible mass, (Mi + M 2 - M)/(M 1 + M 2 ), occurs in the 
equal-mass case and is 1 — l/\/2 « 29%, the same fig- 
ure as obtained by the classical argument involving event 
horizons. In practice, the figure will be much less, since 
M decreases only at the change of topology and is other- 
wise increasing outwards, particularly during the violent 
type-changing phase. 

Here one might prefer on physical grounds to use the 
effective energy E instead of M as a measure of avail- 
able energy, which would require knowledge of angular 
momentum J. One may also evaluate J throughout the 
whole process, but here there is an issue of whether it 
goes smoothly through the topology change, due to the 
more indirect construction and the numerically observed 
spin- flip phenomenon 7 . This is of some interest by it- 
self, as a potential way to relate initial and final angular 
momenta, and could be studied either analytically or by 
pushing simulations inside the region which is normally 
excised. 



XI. CONCLUDING REMARKS 

It seems appropriate to conclude with some important 
recent results and issues. Firstly, there is the uniqueness 
of trapping horizons. Ashtekar & Galloway^ 8 - have shown 
that the structure of a given dynamical horizon is unique, 
in that there is only one way to foliate it by marginal 
surfaces. They also showed that such horizons are not so 
numerous as to foliate a space-time region, so that they 
will tend to interweave one another. 

Andersson, Mars & Simo n 79 ' 80 defined a strictly stable 
marginally outer trapped surface by a condition similar 
to the outer condition and showed that, given such a 



surface in one of a foliation of spatial hypersurfaces, such 
surfaces exist locally in the foliation. This shows that 
trapping horizons are not unique in a given space-time, 
indeed one is locally generated from a given marginal 
surface however the hypersurface is evolved. In practice, 
this is a useful quality since it means that a trapping 
horizon is likely to be found in numerical simulations, 
provided the foliation is chosen reasonably well. 

They also showed that, assuming NEC, the horizon is, 
on any one marginal surface, either spatial everywhere or 
null everywhere. Thus transitions between isolated and 
dynamical phases happen simultaneously with respect to 
the foliation. 

Andersson & Metzger— have shown that, given a spa- 
tial hypersurface with an outer trapped inner boundary 
and an outer untrapped outer boundary, there exists a 
stable marginally outer trapped surface between, which 
is smooth and unique as the outermost such surface. This 
strengthens an earlier result which assumed piecewise 
smoothness^. Since this is the situation expected in 
numerical simulations of black holes, it guarantees the 
existence of an outermost trapping horizon with respect 
to the foliation. 

Another issue is that, in a spherically symmetric space- 
time such as Vaidya, it is possible to find trapped or outer 
trapped surfaces which lie partly outside the spherically 
symmetric trapping horizon, as shown numerically by 
Schnetter & Krishnan^ 3 - and analytically by Ben-Dov^. 
The surfaces are roughly spherical and inside the hori- 
zon except for a long tentacle which extends through it. 
These are rather strange-looking surfaces, so one inter- 
pretation is that, in order to locate a black hole in a 
physically acceptable way, the type of trapping horizon 
needs to be refined further. For instance, it seems rea- 
sonable to require marginal or trapped surfaces to have 
positive Gaussian curvature, which would tend to exclude 
such tentacles. Alternatively or additionally, one might 
require any surface sufficiently close to a marginal surface 
and inside it to be a (non-strictly) trapped surface. 

This review has not addressed quantum issues directly, 
for which one may consult the review of Ashtekar & 
Krishnan 31 . One remarkable feature of the framework 
for isolated horizons is a derivation of black- hole entropy 
proportional to area A. Much less is known about dy- 
namical situations. A very recent result is that, in spher- 
ical symmetry, a tunnelling method derives a local Hawk- 
ing temperature k/2tt precisely for future outer trapping 
horizons^. Thus this framework for dynamical black 
holes seems to be well adapted to the issues raised by 
Hawking radiation. 

Research supported by the National Natural Science 
Foundation of China under grants 10375081, 10473007 
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mission under grant 06DZ111, and by Shanghai Normal 
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